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Abstract. The inclusive and dijet production cross-sections have been measured for jets containing 6-hadrons (6-jets) 
in proton-proton collisions at a centre-of-mass energy of ^J~s = 7TeV, using the ATLAS detector at the LHC. The 
measurements use data corresponding to an integrated luminosity of 34pb~ 1 . The 6-jets are identified using either a 
lifetime-based method, where secondary decay vertices of 6-hadrons in jets are reconstructed using information from 
the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of 
b-hadrons inside jets. The inclusive 6-jet cross-section is measured as a function of transverse momentum in the range 
20 < pt < 400 GeV and rapidity in the range \y\ < 2.1. The 66-dijet cross-section is measured as a function of the 
dijet invariant mass in the range 110 < mjj < 760 GeV, the azimuthal angle difference between the two jets and the 
angular variable x in t wo dijet mass regions. The results are compared with next-to-leading-order QCD predictions. 
Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + 
Pythia. MC@NLO + Herwig shows good agreement with the measured M>-dijet cross-section. However, it does not 
reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta. 



1 Introduction 

The production of fe-quarks in proton-proton collisions at the 
Large Hadron Collider (LHC) provides an important test of 
perturbative QCD (pQCD). Calculations of the fe-quark produc- 
tion cross-section have been performed at next-to-leading order 
(NLO) in pQCD [ ]. These calculations can be matched to dif- 
ferent parton-shower and hadronisation models to produce final 
states that can be compared to those measured in collision data. 

Cross-sections for fo-jet production in high energy pp colli- 
sions have been measured at the SppS [ , ] and Tevatron [4- 
] colliders. The experiments measured cross-sections different 
from those predicted by QCD at the time. This led to substantial 
improvements in the experimental methods and theoretical cal- 
culations. It is therefore of great interest to test the theoretical 
predictions at the higher centre-of-mass energy provided by the 
LHC. Moreover, the measurement of the fe-jet cross-sections is 
an important ingredient in understanding other processes in- 
volving the production of 6-quarks, which represent substantial 
backgrounds in many searches for new physics. Measurements 
of 6-hadron production at s/s = 7TeV in the forward region 
have been reported by LHCb [ ] and in the central region by 
CMS [9, 10]. 

This paper describes measurements of the inclusive 6-jet 
and 66-dijet production cross-sections performed with the AT- 
LAS detector at the LHC. Jets are reconstructed from energy 
clusters in the calorimeter using the anti-fc f algorithm [ ], with 
jet radius parameter R — 0.4. The relatively long lifetime of 
hadrons containing 6-quarks is exploited to obtain a jet sample 



enriched in 6-jets by selecting jets with a reconstructed sec- 
ondary vertex significantly displaced from the primary vertex. 
The number of 6-jets in this enriched sample is derived from 
a fit to the invariant mass distribution of the charged parti- 
cle tracks in the secondary vertex, assuming the pion mass for 
the individual particles. This is referred to as secondary vertex 
mass hereafter. 

The inclusive cross-section is measured for jets containing 
b- or fe-quarks as a function of the transverse momentum, pt, 
and rapidity, y, for jets with 20 < px < 400 GeV and \y\ < 
2.1. The requirement |y| < 2.1 ensures that jets are contained 
within the acceptance of the inner tracking detectors. In the 
kinematic region 30 < p^ < 140 GeV, muon-based 6-tagging 
is used to provide a complementary, and largely independent, 
cross-section measurement as a function of jetpx- 

The 66-dijet cross-section is measured for the leading and 
sub-leading jet in the event as a function of the dijet invariant 
mass, mjj, the azimuthal angle difference between the two jets, 
A(f>, and the angular variable \ = e x P \ui — V2\ for jets with 
Pt > 40 GeV and \y\ < 2.1. The variable x is defined such 
that the cross-section of 2 — > 2 elastic scattering of point-like 
massless particles is approximately constant as a function of 
X — i^cosg* » where 6* is the centre-of-mass scattering angle. 
To measure the cross-sections as a function of \, an additional 
acceptance requirement is used that restricts the boost of the 
dijet system to |yt>oost| = \\v\ + 2/2 1 < 1-1- This reduces the 
sensitivity to parton distribution function (PDF) uncertainties 
at small values of x, where x is the fraction of the proton's mo- 
mentum carried by the parton participating in the hard scatter- 
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ing. The resulting angular distributions provide a test of pQCD 
that is relatively insensitive to PDF uncertainties. 

The measured cross-sections are corrected for all exper- 
imental effects using simulated events, to allow comparison 
with theoretical predictions. 

The data used for these measurements were collected by 
the ATLAS detector in 2010 and correspond to an integrated 
luminosity of 34.0 ± 1.2pb _1 . A detailed description of the 
luminosity determination can be found in Refs. [12, 13]. 



2 The ATLAS detector 

The ATLAS detector [ 1 4] consists of an inner tracking system, 
immersed in a 2 T axial magnetic field, surrounded by electro- 
magnetic calorimeters, hadronic calorimeters and a muon spec- 
trometer. The ATLAS reference system has the origin at the 
nominal interaction point. The x- and y-axes define the trans- 
verse plane, the azimuthal angle (j> is measured around the beam 
axis, z, and the polar angle 9 with respect to the z-axis. The 
pseudorapidity is defined as rj = — In (tan (6/2)). 

The inner detector (ID) has full coverage in <f> and covers 
the pseudorapidity range \<q\ < 2.5. The ID consists of silicon 
pixel and microstrip detectors, surrounded by a transition radia- 
tion tracker (up to |ry| = 2.0). The electromagnetic calorimeter 
is a lead-liquid argon sampling calorimeter covering \rj\ < 3.2. 
Hadronic calorimetry in the barrel (\r)\ < 1.7) is provided by 
a scintillator tile calorimeter using steel as the absorber mate- 
rial. The end-cap hadronic calorimeter uses liquid argon with 
copper absorber plates and extends up to \rj\ = 3.2. Addi- 
tional forward calorimeters extend the calorimetric coverage 
to \rj\ < 4.9, outside the acceptance of this measurement. The 
outer region of the detector is formed by a muon spectrome- 
ter that uses a toroidal magnetic field with a bending power 
of 1.5-5.5 Tm in the barrel and 1.0-7.5 Tm in the end-caps. 
Three layers of muon chambers provide precision tracking in 
the bending plane up to 1 77 1 = 2.7 and the trigger for muons up 
to \r}\ = 2.4. 

The trigger system uses three consecutive trigger levels to 
record a selection of interesting events. The first level trigger 
(LI) is based on custom-built hardware that processes the data 
with a fixed latency of 2.5 us. The second level and the event 
filter, collectively referred to as the high level trigger (HLT), 
are software-based triggers running on computing farms. Then- 
average execution times are 40 ms and 4 s respectively, with a 
design output rate of 3 kHz and 200 Hz respectively. 

Most of the events used in the measurements presented here 
are selected by the calorimeter-based triggers. At LI, the elec- 
tromagnetic and hadronic calorimeters are read out using trig- 
ger towers with a granularity of A<f> x Arj = 0.1 x 0.1, with jet 
identification based on transverse energy in a sliding window 
of 4 x 4 or 8 x 8 trigger towers. At the beginning of data- 
taking in 2010 only the LI triggers were active, while in the 
later runs the HLT was used to refine the jet selection further. 
Events containing jets with 20 < pt < 40 Ge V were triggered 
using the minimum bias trigger scintillators (MBTS) [ ]. The 
MBTS consist of 32 scintillator counters arranged in two discs 
located at ±3.56 m from the interaction point, covering 2.09 < 
\r]\ < 3.84. The hit multiplicity in the MBTS provides a high- 



efficiency trigger for jet events, independent of the jet pr, with 
negligible bias. 

3 Monte Carlo samples and theoretical 
predictions 

Simulated events produced by the Pythia 6.423 [ ] event gen- 
erator are used for the baseline comparisons and to evaluate 
corrections. Pythia implements leading-order (LO) pQCD ma- 
trix elements for 2 —> 2 processes, px-ordered parton-show- 
ers calculated in a leading-logarithmic approximation and an 
underlying event simulation using multi-parton interactions. It 
uses the Lund string model [ ] for hadronisation. All events 
were generated using a specially tuned set of parameters de- 
noted as AMBT1 [ ] with MRST LO* [ ] parton-density 
functions. The generated particles are passed through a full 
simulation [19] of the ATLAS detector and trigger based on 
GEANT4 [ ]. Finally, the simulated events are reconstructed 
and selected using the same analysis chain as is used for the 
collision data, with the same trigger and event selection crite- 
ria. 

The flavour of jets is defined by matching jets to hadrons 
with pt > 5 GeV. The jet is considered a &-jet if a 6-hadron is 
found within AR = yj Atf> 2 + At/ 2 — 0.3 of the jet axis; oth- 
erwise, if a c-hadron is found within the same distance the jet 
is labeled as a c-jet. All other jets are considered light-flavour 
jets. 

The measured cross-sections are compared to NLO pre- 
dictions derived using POWHEG [21-24] and MC@NLO [25, 
26], both using the MSTW 2008 NLO PDFs [2 ] and a 6-quark 
mass of 4.95 GeV. To perform the parton-showering, POW- 
HEG is interfaced to Pythia 6 and MC@NLO to Herwig 6 [28]. 
For Herwig, the AUET1 [ ] tune is used. In contrast to Pythia, 
Herwig uses an angular-ordered parton-shower model and a 
cluster hadronisation model. 



4 Event and jet selection 

The events used in the lifetime-based analysis are triggered by 
the LI or HLT jet triggers, with the exception of the 20 < 
Pt < 40 GeV bin in the inclusive cross-section measurement 
where the MBTS trigger is used. The trigger efficiency for b- 
jets using these trigger selections is estimated to be above 97% 
in all cases and typically close to 100%. For the muon-based 
cross-section measurement the combination of a jet and a muon 
trigger is required, which results in an efficiency ranging from 
about 35% for jets with pt < 50 GeV to 65% for jets with 
Pt > 105 GeV. While this efficiency is lower, the different 
trigger prescale factors allocated a much higher rate to the jet- 
muon trigger than to the inclusive jet triggers for a similar jet 
Pt threshold. 

Quality selections are applied to the reconstructed jets to 
ensure that they are not produced by poorly calibrated detec- 
tor regions or noisy calorimeter cells [30]. Additionally, the 
charged particle tracks contained in the jets are required to be 
of adequate quality for ^-tagging [ ] and a good reconstructed 
primary vertex is required that contains at least 10 tracks with 
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Fig. 1. Examples of template fits to the measured distribution, before and after applying the requirement of L/ol > 5.85. The error bars 
represent the data statistical errors. The differences between the data and the sum of the templates are covered by the systematic uncertainties 
on the template shapes. 



Pt > 150 Me V. The combined efficiency of the reconstruction 
and the quality requirements is determined to be above 96% for 
6-jets. 

The secondary vertex 6-tagging algorithm used, SV0 [31], 
aims at reconstructing the position of the displaced vertex from 
the charged decay products of long-lived particles in a jet. The 
SV0 algorithm reconstructs two-track vertices from tracks in- 
side a cone of AR = 0.4 around the jet axis that are signifi- 
cantly displaced from the primary vertex, based on the three-di- 
mensional impact parameter significance. Quality requirements 
are applied to the two-track vertices to reject vertices that are 
compatible with the primary vertex, are located at a radius con- 
sistent with one of the pixel detector layers or contain tracks 
that have an invariant mass consistent with a K s meson, a A 
baryon or a photon conversion. A single secondary vertex is 
then fitted to all the tracks which contribute to any of the re- 
maining two-track vertices in the jet. 

The signed decay length significance of the secondary ver- 
tex, L/ctl, is used to select a jet sample enriched in 6-jets. 
The sign of the decay length is given by the sign of the pro- 
jection of the decay length vector onto the jet axis. Jets with 
L/<jl > 5.85 are referred to as 6-tagged jets. The selection at 
5.85 is chosen such that it produces a 50% 6-tagging efficiency 
for 6-jets in simulated ti events. 



4.1 6-tagging efficiency 

The efficiency of the chosen selection on L/<jl is estimated 
with a data-driven method that uses jets containing a muon. The 
number of 6-jets before and after 6-tagging can be obtained us- 
ing the variable , which is defined as the momentum of the 
muon transverse to the combined muon plus jet axis. Muons 
originating from 6-hadron decays have a harder p 1 ^ 1 spectrum 
than muons in c- and light-flavour jets. Templates of the p^: 1 
shape are constructed for each jet flavour separately. The tem- 



plates for 6- and c-jets are extracted from Monte Carlo simula- 
tion, while the light-flavour template is obtained from a light- 
jet enriched data sample. These are then fitted [ ] to the mea- 
sured p 1 ^ 1 spectrum of muons in jets to obtain the fraction of 
6-jets before and after requiring a 6-tag. The fit determines the 
relative contributions of the 6-, c- and light-flavour templates 
such that their sum best describes the shape of the p 1 ^ 1 distri- 
bution in data. Having obtained the flavour composition of jets 
containing muons from the pJj? 1 fits, the 6-tagging efficiency is 
defined as 

j-tag _ Artag 
c data _ Jb r> 

where and /^ ag are the fractions of 6-jets before and after 6- 
tagging is applied, and N and A ftag are the total number of jets 
in those two samples. The factor C corrects the efficiency for 
biases introduced by differences between data and simulation 
in the modelling of the 6-hadron direction and by heavy -flavour 
contamination of the p r ^ 1 template for light-flavour jets. The 
magnitude of these corrections is typically a few percent. Ex- 
amples of fits to the p 1 ^ 1 distribution before and after the L/a^ 
requirement are shown in Fig. 1. 

The p™ 1 method can be used to determine the 6-tagging ef- 
ficiency for 6-jets containing 6-hadrons that decay semileptoni- 
cally. Studies have been performed to show that this determina- 
tion can be extended to all 6-jets and a systematic uncertainty 
due to this generalization is assigned to the 6-tagging efficiency 
for all 6-jets. A detailed account of the systematic uncertainties 
in the 6-tagging efficiency calibration is given in Ref. [33]. 

The discriminating power of the p 1 ^ 1 method decreases with 
increasing jet px> hence this method can only provide a data- 
driven determination of the 6-jet tagging efficiency for jet px 
values up to about 140 GeV. For jets with px > 140 GeV, the 
6-tagging efficiency is derived from simulation and multiplied 
by a correction factor of 0.88 ± 0.18 that accounts for the dif- 
ference between data and simulation observed in the pt range 
90-140 GeV. Comparisons between data and simulation as a 
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Fig. 2. Examples of purity fits in the inclusive and dijet measurements. The error shown for the 6-fraction is the uncertainty on the fit parameter. 
For the inclusive measurement the statistical uncertainty on the sum of the templates, indicated by the shaded area, is taken into account in 
the fit. In the dijet measurement the templates are parameterized, the uncertainty on the parameterization is taken into account as a systematic 
uncertainty and not shown here. 



function of jet px show that the simulation models the data 
equally well in all regions of jet px in which data measure- 
ments are available, so the above extrapolation is well moti- 
vated. Moreover, detailed comparisons between data and sim- 
ulation as a function of jet px, in terms of the quantities that 
affect 6-tagging, show that the effect of any mismodelling of 
the 6-tagging performance at higher jet px values is within the 
systematic uncertainties assigned to the 6-tagging efficiency. 
The efficiency after applying the requirement of L/a^ > 5.85 
ranges from 20% for 6-jets of px < 40 GeV and \y\ > 1.2 to 
55% for central 6-jets with px of about 100 GeV. 

The p 1 ^ 1 distribution can also be used as a discriminant vari- 
able to measure the inclusive fo-jet cross-section directly. While 
this method is statistically limited and cannot be used beyond 
140 GeV, as mentioned above, it does provide a useful cross- 
check for the lifetime-based measurement. Many of the sys- 
tematic uncertainties are different and the sample of jets used 
is statistically largely independent from that used in the life- 
time-based measurement. The muon-based cross-section mea- 
surement is described in Section 5. 



4.2 £>-jet purity 

In the lifetime-based measurement, the fraction of 6-jets in the 
6-tagged sample of jets, referred to as the purity of the sample, 
is determined by performing a template fit to the secondary ver- 
tex mass distribution. The templates for b-, c- and light-flavour 
jets are extracted from Monte Carlo simulation. The average in- 
variant mass of a secondary vertex increases when going from 
light-flavour jets via c-jets to 6-jets, making it possible to sep- 
arate the flavours by determining the relative fractions of the 
templates that best describe the vertex mass distribution in data. 

For the inclusive cross-section measurement, the number 
of b-, c- and light-flavour jets is fitted by maximizing a binned 



likelihood function that takes into account the statistical uncer- 
tainties in both the data and the templates. The fit is performed 
for each px and y region separately, in vertex mass bins of 
200 Me V. 

In the dijet cross-section measurement, the fraction of 6-jet 
pairs is determined from a template fit to the sum of the vertex 
masses of the two 6-tagged jets. This fit uses two templates: 
the 6-template, where both jets are matched to a 6-hadron in 
simulation; and a non-6 template, where at least one of the two 
jets is a c- or light-flavour jet. In order to reduce the effect of 
the limited statistics in simulation, a parameterization is used to 
smooth the templates. The fit is performed for each kinematic 
region separately. Typical fit results in the inclusive and dijet 
measurements are shown in Fig. 2. 

5 Results 

All the measured cross-sections are corrected for experimental 
effects using a bin-by-bin correction, so as to represent parti- 
cle-level cross-sections of jets containing 6-hadrons. The cor- 
rection is obtained from Pythia simulated dijet events by calcu- 
lating the cross-sections for both particle-level 6-jets (including 
muons and neutrinos) and reconstructed 6-jets. The correction 
factors are derived bin-by-bin in each distribution by taking the 
ratio of the two cross-sections. 



5.1 Systematic uncertainties 

The dominant systematic uncertainties in both the inclusive and 
the dijet cross-section measurements come from the fo-jet en- 
ergy scale calibration, and the determination of the 6-tagging 
efficiency and purity. The systematic uncertainties, including 
those on the muon-based measurement which will be discussed 
in Section 5.2, are summarized in Table 1. 
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Table 1. Summary of the most important systematic uncertainties on 
the lifetime-based inclusive 6-jet and 66-dijet cross-section, and on the 
muon-based cross-section measurement. 



Syst. uncertainty 


Inclusive 6-jet 


66-dijet 


Muon-based 


Jet energy scale 


10-20% 


10-20% 


15-20% 


6-tagging efficiency 


5-20% 


30-50% 




6-jet purity fit 


3-8% 


20-30% 


8-18% 


Luminosity 


3.4% 


3.4% 


3.4% 


Other sources 


2% 


2% 


3% 



Jets are calibrated to the hadronic scale using the inclusive 
jet energy scale calibration [34, 35], which is based on p?- and 
^-dependent correction factors derived from Monte Carlo sim- 
ulation and validated with test beam measurements. The un- 
certainty on this jet energy scale varies between 2% and 6% 
depending on the jet px and rapidity region. 

For heavy-flavour jets, two studies were performed to es- 
timate additional contributions to the jet energy scale uncer- 
tainty that account for flavour-dependent systematic uncertain- 
ties. Firstly, the uncertainty on the calorimeter response for 6- 
jets due to their different particle composition has been evalu- 
ated using single hadron response studies [ ]. This method 
compares the relative response of 6-tagged jets in ti events 
with that of inclusive jets in QCD dijet events. For jets within 
\i]\ < 0.8 and 20 < px < 250 GeV, this difference is found to 
be negligible (<0.5%). Secondly, systematic uncertainties for 
6-jets were studied in Monte Carlo simulation by comparing 
particle-level jets to reconstructed jets. The variations that were 
studied include the modelling of fragmentation, hadronisation, 
parton-showers and the underlying event, but also variations in 
soft-physics tunes and the effects of the uncertainty on the ma- 
terial description. The 6-jet energy scale uncertainty obtained 
using these two methods is validated in data by comparing the 
total transverse momentum of the calorimeter jet to that of the 
charged particle tracks associated to it [37]. 

It is found that there is an additional 2.5% uncertainty on 
the 6-jet energy scale with respect to the uncertainty on the 
energy scale of inclusive jets. This extra uncertainty is added 
in quadrature. When propagated to the cross-section measure- 
ments, this leads to an uncertainty of 10% to 20%, depending 
on the kinematic region. 

The most important contributions to the systematic uncer- 
tainty on the 6-tagging efficiency originate from the modelling 
of muons in jets in the simulation, the generalization of the 
efficiency from 6-jets with muons to inclusive 6-jets, and the 
limited statistics of the templates used for the p^ 1 fi ts - More 
details about the 6-tagging efficiency uncertainty can be found 
in Ref. [33]. The resulting uncertainty on the cross-sections 
amounts to between 5% and 20% for the inclusive 6-jet cross- 
section, and between 30% and 50% for the dijet cross-section. 

The systematic uncertainties from the purity fits account 
for the observed differences between jets in collision data and 
those in the Monte Carlo simulation used to derive the tem- 
plates. The uncertainty is derived from studies of the secondary 
vertex mass distribution in light-jet enriched samples and 6-jet 
enriched samples. The light-jet enriched sample is obtained by 
selecting jets with a negative decay length. For the 6-jet en- 



riched samples, two methods are used: the first requires another 
6-tagged jet to be present in the event, while the second selects 
secondary vertices with high track multiplicities. The observed 
differences in the secondary vertex mass distribution are then 
used to correct the template shapes and re-evaluate the fits. The 
difference in the cross-section is found to be between 3% and 
8% and this is assigned as the purity fit systematic uncertainty. 
For the 66-dijet cross-section, the most important contribution 
to the systematic uncertainty on the 66-fraction is due to the 
limited template statistics. The effect of the statistical uncer- 
tainty of the templates is estimated by varying the shape param- 
eters of the parameterized templates within their uncertainties 
and re-evaluating the cross-section. The resulting uncertainty 
is between 20% and 30%. 

The systematic uncertainty on the luminosity determina- 
tion is 3.4% [ ]. The remaining sources of systematic uncer- 
tainty, such as the effect of possible differences in the cross- 
section shapes between data and simulation on the bin-by-bin 
corrections, differences in the jet energy resolution between 
data and simulation, the trigger efficiency and the jet selection 
efficiency, lead to a combined systematic uncertainty of about 
2%. 

The effect of different shower and hadronisation models 
is included in the jet energy scale uncertainty. The impact of 
changing the shape of the px distribution on the bin-by-bin 
corrections was found to be much less than 1%. Using Herwig 
instead of Pythia to derive the correction factors gives statisti- 
cally consistent results. 

5.2 Muon-based 6-jet cross-section 

The p!^ 1 method, used for calibrating the 6-tagging efficiency, 
is also used to obtain an independent measurement of the in- 
clusive 6-jet cross-section in the range 30 < px < 140 GeV. 
This measurement uses jets containing a muon of px > 4 GeV 
within a cone of AR — 0.4 from the jet axis. The flavour 
composition of this sample is extracted from a template fit to 
the muon p^ 1 distribution. The templates for 6- and c-jets are 
obtained from Monte Carlo simulation. Two data-driven tech- 
niques are employed to extract the shape of the muon p!^ 1 in 
light-flavour jets. The first takes the shape from jets with neg- 
ative decay length in data, which is then corrected for 6-jet 
contamination using simulation. The second method uses in- 
clusive jets without a muon; the template is then obtained by 
converting each track inside the jet into a muon and weighting 
the resulting pJj? 1 by a probability to simulate hadron decays in 
flight. The 6-jet fraction is evaluated using both methods, tak- 
ing the average as the central value and assigning the difference 
between them as a systematic uncertainty. 

The dominant sources of systematic uncertainties in this 
measurement are the 6-jet energy scale (15-20%) and the purity 
fits (8-18%). Contributions to the purity fit systematics include 
limited template statistics and uncertainties in the modelling 
of semileptonic 6-hadron decays and 6-fragmentation. The first 
modelling error is estimated by varying the muon momentum 
distribution in the rest frame of the 6-hadron between that mea- 
sured by DELPHI [38] and that measured by BABAR [39]. 
The second is measured by varying the fraction of the 6-jet en- 
ergy carried by the 6-hadron by ±5% and rederiving the 6-jet 
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templates in the simulation. Apart from the 6-jet energy scale, 
the systematic uncertainties are to a large extent specific to the 
muon-based measurement. This makes the comparison with 
the lifetime-based cross-section measurement a useful cross- 
check. 



5.3 Cross-section results and discussion 

The double-differential inclusive 6-jet cross-section is shown in 
Fig. 3 as a function of jet in four different rapidity regions. 
Figure 4 shows the single differential cross-section as a func- 
tion of px> integrated over the entire rapidity range of \y\ < 
2.1. In the px range where the lifetime-based and the muon- 
based measurements overlap, both results are shown. The data 
are compared to NLO predictions derived with POWHEG and 
MC@NLO. In addition, the data are compared to the Pythia 
prediction. Pythia, as a leading-logarithmic parton-shower gen- 
erator, is not expected to predict the correct normalization. The 
Pythia prediction is scaled by a factor x0.67 in order to match 
the measured integrated cross-section, allowing a comparison 
of the cross-section shapes. All three calculations describe the 
general features of the cross-section reasonably well. 
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Fig. 3. Inclusive double-differential 6-jet cross-section as a function 
of pt for the different rapidity ranges. The data are compared to the 
predictions of Pythia, POWHEG and MC@NLO. The leading-order 
Pythia prediction is scaled (x0.67) to the measured integrated cross- 
section. 



To allow for a better comparison between the data and the 
NLO predictions, Fig. 5 shows the ratio of the measured cross- 
section to the NLO theory predictions for \y\ < 2.1 (top) and 
for each rapidity region separately. The plot for the full rapidity 
acceptance also allows a direct comparison between the life- 
time-based and the muon-based cross-section measurements 
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Fig. 4. Differential 6-jet cross-section as a function of pt for 6-jets 
with \y\ < 2.1. The data are compared to the predictions of Pythia, 
POWHEG and MC@NLO. In the region 30 < p T < 140 GeV the 
muon-based cross-section measurement is also shown. For the muon- 
based measurement only the POWHEG prediction is shown. 



in the overlapping range, indicating a good agreement be- 
tween the two measurements. Good agreement is also observed 
between the measured cross-section and the NLO predictions 
obtained using POWHEG + Pythia in all rapidity regions. MC- 
@NLO + Herwig, however, predicts a significantly different 
behaviour of the double-differential cross section, as shown in 
Fig. 5b. When the cross-section is integrated over the full ra- 
pidity acceptance this effect averages out somewhat and MC- 
@NLO + Herwig shows better agreement with data. It has been 
checked that the qualitative behaviour remains the same when 
POWHEG is interfaced to Herwig instead of Pythia, implying 
that the observed rapidity dependence in MC@NLO + Herwig 
is not resulting from the parton-shower Monte Carlo program. 
On the other hand, POWHEG + Herwig appears to predict a 
cross-section that is consistently lower than the POWHEG + 
Pythia prediction. This would suggest that the deficit of MC- 
@NLO + Herwig compared to the data in Fig. 5b, may be 
partly due to the Herwig parton-showering. 

Comparison to the inclusive (all-flavour) jet cross-section 
measurement [ ], shows that the fraction of jets containing a 
6-hadron is approximately 5% in the kinematic region where 
the two measurements overlap, 60 < j»t < 400 GeV and \y\ < 
2.1. 

The 66-dijet cross-section is shown as a function of dijet 
mass in Fig. 6. It should be noted that nearby 66-pairs, as ex- 
pected for example from gluon splitting, are generally not re- 
solved as separate jets. Also, since the measurement refers to 
the leading and sub-leading jet in the event, the contribution 
from gluon splitting is expected to be small. The bfe-dijet cross- 
section is compared to Pythia and the NLO predictions ob- 
tained using POWHEG and MC@NLO. The Pythia prediction 
is again normalized to the measured integrated cross-section, 
here using a factor of x0.85. The Pythia normalization is not 
expected to be the same as that used in the inclusive cross- 
section, given the different event selection used. All theory pre- 
dictions show good agreement with the measured cross-section. 

Figure 7 shows the fractional fe5-dijet cross-section as a 
function of the azimuthal angle between the two jets, A<p. The 
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measurement overlaps with the muon p%r measurement both results are shown. The top plot shows the full rapidity acceptance, while the four 
smaller plots show the comparison for each of the rapidity ranges separately. The data points show both the statistical uncertainty (dark colour) 
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the statistical uncertainty only. Systematic uncertainties in the NLO predictions are discussed in the text. 



dijets selected in this measurement show a pronounced back- 
to-back: configuration in the transverse plane that is generally 
well reproduced by QCD generators. 

The fefe-dijet cross-section as a function of the angular vari- 
able x 1S shown in Fig. 8 for dijets with |yboost| < 1.1. The \ 
distribution is well reproduced by the theoretical calculations. 
The distribution flattens for large invariant mass values. 

In the NLO calculations, the renormalization and factoriza- 
tion scales are set equal to the transverse energy of the hard- 



est parton: Q~ 



p^. To estimate the potential 



impact of higher order terms not included in the NLO calcu- 
lation on the theory predictions, the renormalization scale is 
varied from half to twice its default value. Similarly, to esti- 
mate the impact of the choice of the scale where the PDF evo- 
lution is separated from the matrix element, the factorization 
scale is varied up and down by a factor of two. The effect of 
each of these variations on the NLO cross-section prediction 
is estimated using POWHEG and found to be approximately 
20% for all kinematic regions. Finally, the uncertainty on the 
PDFs is estimated by deriving the NLO predictions using the 
NNPDF [40] and CTEQ 6.6 [4 ] PDFs, resulting in a differ- 
ence of approximately 10% for all kinematic regions. 



6 Conclusions 

The inclusive 6-jet and 66-dijet production cross-sections have 
been measured in proton-proton collisions at a centre-of-mass 
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energy of 7TeV, using data with an integrated luminosity of 
34 pb -1 recorded by the ATLAS detector. 

The inclusive fo-jet cross-section was measured as a func- 
tion of jet px in the range 20 < pt < 400 GeV and rapidity in 
the range \y\ < 2.1. The fofo-dijet cross-section was measured 
as a function of dijet invariant mass in the range 110 < mjj < 
760 GeV, as a function of the azimuthal angle difference and 
of the angular variable \- The measurements are dominated 
by systematic uncertainties, mainly coming from the fo-jet en- 
ergy scale and the determination of the fo-jet tagging efficiency 
and purity. The measured cross-sections have been compared 
to next-to-leading order QCD predictions derived using POW- 
HEG interfaced to Pythia and MC@NLO interfaced to Herwig. 

The inclusive cross-section measured over \y\ < 2.1 for 
fo-jets identified by the presence of a secondary vertex is com- 
pared to a largely independent cross-section measurement that 
uses muon-based fo-tagging in the range 30 < pr < 140 GeV. 
The two measurements show good agreement. 

The inclusive fo-jet cross-section is found to be in good 
agreement with the POWHEG + Pythia prediction over the full 
kinematic range. MC@NLO + Herwig, however, predicts a sig- 
nificantly different behaviour of the double-differential cross 
section that is not observed in the data. The normalized leading- 
order Pythia prediction shows broad agreement with the mea- 
sured cross-section. 

POWHEG + Pythia and MC_@NLO + Herwig show good 
agreement with the measured fofo-dijet cross-sections, as does 
the normalized leading-order Pythia generator. 
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